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Abstract: A phenomenological analysis of radiative V — * P7 and P — * V"f decays is per- 
formed with the purpose of determining the gluonic content of the r/ and 77' wave functions. 
Our results show that within our model there is no evidence for a gluonium contribution in 
the rj, = 0.00 ±0.12, or the 77', Z£ = 0.04±0.09. In terms of a mixing angle description 
this corresponds to cfip = (41.4±1.3)° and l^n'Gl = (12 ±13)°. In addition, the rj-rj' mixing 
angle is found to be (ftp = (41.5 ± 1.2)° if we don't allow for a gluonium component. 
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1. Introduction 

Very recently the KLOE Collaboration has reported a new measurement of the ratio 
= B{4> —* ffl)/B{<f> — * 777) 0]- Combining the value of R<f, with other constraints, 
they have estimated the gluonium content of the rj meson as = 0.14 ± 0.04, which 
points to a significant gluonium fraction in the rf wave function incompatible with zero by 
more than 3<7. This new result contrasts with the former value Zj = 0.06±g;g§, which was 
compatible with zero within la and consistent with a gluonium fraction below 15% 0|. 
Both experimental analysis obtain their results using the same set of constraints derived 
from other measured ratios, namely 1 Y(r( — > 77)/T(7r° — > 77), T(r]' — > p^)/T(uj — > tt°j) 
and r(rj' — > oj~j)/T{oj — > vr°7), as well as R^. The sole difference between the two analyses, 
apart from the obvious improvement in the precision of the new measurements, is the 
inclusion in the amplitudes of Ref. |]|] of two extra parameters to deal with the overlap of 
the vector and pseudoscalar meson wave functions produced in the transitions V — * P7 or 
P — > 1/7, a feature first introduced in Ref. Q. However, the new analysis of Ref. Q uses 
the most recent experimental data taken from Ref. H in association with the values for 
the parameters related to the overlap which were obtained in Ref. Q from a fit to available 
experimental data at that time. Therefore, a reanalysis of this uncomfortable situation is 
needed before drawing definite conclusions on the gluon content of the r/ and rf mesons. 
The purpose of this work is to perform a phenomenological analysis of radiative V — > P7 
and P — > V^f decays, with V = p,K*,u;,(j) and P = n,K,r},rf, aimed at determining 

lr rhe analysis in Ref. B did not contain the constraint imposed by V(f]' — ► lo^)/T{lo — » tv°j). Neverthe- 
less, the inclusion of this constraint does not change the main results of the analysis, as seen in Ref. 



- 1 - 



the gluonic content of the 77 and rf wave functions. Similar analysis were driven in the 
seminal work by Rosner H, where the allowed gluonic admixture in the r\ was found to 
be small, \Z^\ < 0.4, but not necessarily zero. The gluonic content of the rj' could not be 
established due to the lack of data on (ft — > 7/7. Later on, Kou pointed out that the 7/ 
gluonic component might be as large as 26% ||. These two works also included in their 
analysis the constraints provided by P —* 77 with P = rj,r]'. Here, we will not take into 
account these decays since they are nicely described in a two mixing angle scheme, as shown 
in Refs. |@, §■ A more recent analysis studying the pseudoscalar glueball-qq mixing from 
J/V#') -> VP decays gets \Z V \ = 0.042 and \Z V >\ = 0.161 §. 

In Sec. E|, we briefly introduce the notation. The model for VPj Ml transitions to 
be used in our analysis is presented in Sec. ||[ This model includes the overlap of vector 
and pseudoscalar wave functions as a main ingredient |J. As stated above, the use of old 
values for the overlap parameters in addition to the latest data may be at the origin of the 
discrepancy concerning the gluon content of the rf between the two KLOE analyses [l], |2|. 
Sec. |4] is devoted to the data fitting of the most recent V — ► P7 and P — ► V7 experimental 
data with the aim of finding the gluonic admixture in the 77 and rf wave functions. As a 
matter of comparison with other approaches, a graphical representation of our main results 
is shown in Sec. |. Finally, in Sec. ^ we summarize our work and present the conclusions. 
A decomposition of the mixing parameters in terms of Euler angles is left for App. |A]. 

2. Notation 

We will work in a basis consisting of the states \rj q ) = ~^\ U ^ L + dd), \rj s ) = \ss) and 
\G) = |gluonium). The physical states rj and rf are assumed to be linear combinations of 
these 

\ V )=X v \r jq )+Y v \ Vs ) + Z v \G) , 
\rf) = X v ,\r] q ) + Y v ,\rj s ) + Z V ,\G) , 
with Xl ( Zl ( = 1 and thus Xl ( +Yl t) < 1. A significant gluonic admixture 



in a state is possible only if Z\ ,-, = 1 — X\ ,-, — Y^< > [^]. This mixing scheme 



v(v') viv') v(v') v(v') nW) 

assumes isospin symmetry, i.e. no mixing with ir , and neglects other possible admixtures 
from cc states and/or radial excitations. In absence of gluonium, Z^i\ = 0, the mixing 
parametrization fl2,l| ) is reduced to the standard pattern in the quark-flavour basis 

1 77) = cos^pI??,) -sin^p|7? s ) , 

(2.2) 

\t]'} = s'm<p P \ri g } +cos<pp\r] s } , 

with X v = Yff = cos 4>p, Xrf = -Y v = sin^p, and X*,,s + Y*,,-, = 1. Similarly, for the 
vector states u and <f> the mixing is given by 

\uj) = cos 4>v\u q ) — sm4>v\4> s ) , 

(2.3) 

\4>) = smc/) V \u; q ) + cos (j) V \(j) s ) , 

where \cv q ) and \4> s ) are the analog non-strange and strange states of \rj q ) and |?7 S ), respec- 
tively. 
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3. A model for VPj Ml transitions 



We will work in a conventional quark model context and assume that pseudoscalar and 
vector mesons are simple quark-antiquark S"-wave bound states. All these hadrons are thus 
extended objects with characteristic spatial extensions fixed by their respective quark- 
antiquark P or V wave functions. In the pseudoscalar nonet, P = TT,K,r],r]', the quark 
spins are antiparallel and the mixing pattern is given by Eq. In the vector case, 

V = p, K* , u, <j), the spins are parallel and mixing is similarly given by Eq. (|2.3| ). We will 
work in the good SU{2) limit with m u = = fh and with identical spatial extension of 
wave functions within each P and each V isomultiplet. SU(3) will be broken in the usual 
manner taking constituent quark masses with m s > fh but also, and this is a specific feature 
of our approach, allowing for different spatial extensions for each P and V isomultiplet. 
Finally, we will consider that even if gluon annihilation channels may induce rj-rj' mixing, 
they play a negligible role in VP'y transitions and thus fully respect the usual OZI-rulc. 

In our specific case of VP'y Ml transitions, these generic statements translate into 
three characteristic ingredients of the model: 



i) A VP'y magnetic dipole transition proceeds via quark or antiquark spin- flip ampli- 
tudes proportional to p q = e q /2m q . Apart from the obvious quark charge values, 
this effective magnetic moment breaks SU(3) in a well defined way and distinguishes 
photon emission from strange or non-strange quarks via m s > fh. 



ii) The spin-flip V <-> P conversion amplitude has then to be corrected by the relative 
overlap between the P and V wave functions. In older papers [1C, [ll|] a common, 
flavour-independent overlap was introduced. Today, with a wider set of improved 
data, this new symmetry breaking mechanism can be introduced without enlarging 
excessively the number of free parameters. 



Hi) Indeed, the OZI-rule reduces considerably the possible transitions and their respective 
VP wave-function overlaps: C s , C q and C w characterize the (r] s \(j) s ), (r] q \Lo q ) = {rj q \p) 
and {i^\ijJ q ) = spatial overlaps, respectively. Notice that distinction is made 

between the \tt) and \rj q ) spatial extension due to the gluon or U (1)a anomaly affecting 
the second state. Independently, we will also need Ck for the (K\K*) overlap between 
strange isodoublets. 

It is then a trivial task to write all the VP'y couplings in terms of an effective g = gui q -wY- 

g p o^o 1 = g P +7r+7 = \g , g^ = g cos (p v , g^ = g sin <f) V , 

9k*o k ^ = -\gz K (l + ^) , 9k*+k+ 1 = \gZK (2- ^7) , (3.1) 

Qprj^ = 9 z q A?} , 9pri'~j = 9 %q -^V > 



- 3 - 



dujiry = \g [z q X v cos <p v + 2^z s Y n sin (f> v 
3w»j'7 = §5 (z q x v > cos </>y + 2^z s Yrf sin </>y 



= ^9 { z q X v sax <j) V - 2^-z s Y v cos <f) V 



9<Wt = \s {z q X v > sin </>y - 2^-z s Y v , cos 0y) , 

where we have redefined z q = C q /C n , z s = C s jC v and zk = Ck/C^. The normalization 
of the couplings is such that g W7r y = g cos <pv = 2 + /jj) cos (f>v = eC,,- cos (fry /m and 
the decay widths are given by 



^-P7)44?IP/. r(P^ 7 ) = 4?|p 7 |*, (3.3) 



where p 7 is the final photon momentum. 
4. Data fitting 

We proceed to fit our theoretical expressions for the amplitudes in Eqs. ( |3.1| , |3.2D comparing 
the available experimental information on r(V — > P7) and r(P — > V7) taken exclusively 
from Ref. [|| with the corresponding decay widths in Eq. ( |3.3| ). Looking at these amplitudes, 
one immediately realizes that the overlapping parameters z q>s and the mixing parameters 
X v ( v i\ and Y v t v ^ always appear in pairs, namely, ZqX^M) and z s Y v ^ v iy Constraining 
these four different combinations will fix four parameters at most. However, there are five 
independent parameters to be fixed from them, viz. z q)S and the three mixing parameters 
related to the most general case of accepting a gluonic admixture in both the 77 and rf 
mesons (see below). So, either we fix the z's to unity and then constrain the three mixing 
parameters of the general case, or we leave the z's free and then we are restricted to allow 
for gluonium either in the rj or 7/ wave function only. In the following, we consider these 
possibilities. Both are interesting since a comparison of their results will allow us to check 
the relevance of taking into consideration the overlapping parameters, which are specific of 
our approach. Furthermore, leaving the z's free will permit us to fix the gluonic content 
of the rf in a way identical to the experimental measurement by KLOE, that is, under the 
hypothesis of no gluonium in the rj wave function. Unfortunately, due to the pairing of 
parameters mentioned above, a simultaneous fit of the z's and the gluonic admixture in 
the rj and rf is not possible. 

As shown in detail in App. |A[ due to the orthonormality conditions in Eqs. ( A. 2 , A.3| ) 



the mixing pattern of rj and rf is described by means of three independent parameters, 
which we choose to be X v , X rj i and Y r] i . Accordingly, the parameters related to the gluonic 
content, and Z„/, and Y^ are written in terms of the former ones as 



\ Z v,(v')\ — V 1 X vM) Y vM) 



X v X Tf Y r ,, + J(l - X* - Y*)(l - XI - Xj) 



(4.1) 



^ 1 — X'^i 
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The mixing parameters can also be expressed in terms of three angles, 4>p, ip^G and 4> v 'g, 
the two latter weighting the gluonic admixture in the rj and rj , respectively. Using this 
angular parametrization one gets 

X„ = cos 4>p cos (j)^ , Xrf = sin <f>p cos </y G — cos <f>p sin sin ^g , 

Y r] = — sin <f>p cos 4>rjG , Yrf = cos 4>P cos <?Vg + sm 4>P Sin 4>nG sin 4>^g > (4-2) 

Z v = - sin ^ G , Z v > = - sin (/y G cos ^ . 

Fits to experimental data are performed imposing the constraints in Eq. ( |4.1| ) or, equiva- 



lently, using the decomposition in Eq. (4.2) 



We start considering the first of the possibilities noted before. Thus, we assume that 
the overlap of the P and V wave functions is flavour-independent, i.e. C q = C s = Ck = CV 
and hence z q = z s = z% = 1. The fit in this case is very poor, x 2 /d.o.f.=31.2/6. The 
quality of the fit gets worse when (J) v g and 4> n 'G are set to zero, x 2 /d.o.f.=45.9/8 with 
<f>P = (41.1 ± 1.1)°. 

Clearly, in order to obtain a good fit one has to relax the constraint imposed on the 
overlapping parameters. Hence, we begin to discuss the second of the possibilities, that is 
to say, to leave the z's free and restrict the gluon content of the rj or rj meson. However, 
as a matter of comparison, we first consider the absence of gluonium in both mesons, 
i.e. 4> v g = (ftr/'G = 0- I 11 addition, we also fix the vector mixing angle <j)y to its measured 
value tan cj) V = +0.059 ± 0.004 or cj) V = (3.4 ± 0.2)° p| and the ratio of constituent quark 
masses to m/m s ~ 1/1.45. The fit in this case is not yet satisfactory, x 2 /d.o.f.=14.0/7. 
The quality of the fit improves when the ratio m/m s is left free, x 2 /d.o.f.=7.6/6 with 
m s /fh = 1.24 ± 0.07. If 4>y is also left free, the final result of the fit gives x 2 /d.o.f.=4.4/5 
with 

g = 0.72 ± 0.01 GeV" 1 , <f> P = (41.5 ± 1.2)° , <j) V = (3.2 ± 0.1)° , 

(4.3) 

1.24 ±0.07, z q = 0.86 ± 0.03 , z s = 0.78 ±0.05, z K = 0.89 ± 0.03 . 



in 



The fitted values for the two mixing angles <f> P and <t>y are in good agreement with most 
results coming from other analyses using complementary information (see, for instance, 



Ref. 1 13] and references therein). Our value for the pseudoscalar mixing angle also agrees 



with the latest measurement from KLOE, <pp = (41.4 ± 1.0)° Q. The free parameters z's 



are specific of our approach and are not fixed to one as in previous analyses 1C, 11]. As 
mentioned, if we fix the z's to unity, the fit gets much worse (x 2 /d.o.f.=45.9/8). This shows 
that allowing for different overlaps of quark-antiquark wave functions and, in particular, 
for those coming from the gluon anomaly affecting only the rj and rj singlet component, is 
indeed relevant. 

Finally, the symmetry breaking parameter m s /rh which is essential to adjust the ratio 
between the two K*-K transitions is affected by large uncertainties due to the difficulty 
in extracting the neutral and charged K* — > K^y widths from Primakoff-effect analyses 2 . 



2 The neutral and charged K* —* Kj transitions have been measured only by one and two experimental 
groups, respectively, and seem to need further confirmations M. 
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Transition 


5 ^ 7 (PDG) 


gf Pl {Yit 1) 


3&, 7 (Fit 2) 


P ~ 


-> m 


0.475 ± 0.024 


0.461 ±0.019 


0.464 ± 0.030 


rf- 


♦ A 


n a i 4- n nQ 


n a i 4_ n no 
U.41 ± u.uz 


n ac\ 4- n n/i 

U.4U ± U.U4 


UJ - 


-> 777 


0.140 ±0.007 


0.142 ±0.007 


0.143 ±0.010 


rf- 


-> 


0.139 ±0.015 


0.149 ±0.006 


0.146 ±0.014 


</>- 


-> r/7 


0.209 ± 0.002 


0.209 ± 0.018 


0.209 ±0.013 


4>- 


* »/7 


0.22 ±0.01 


0.22 ±0.02 


0.22 ±0.02 



Table 1: Comparison between the experimental values <?™p 7 (in GeV x ) for the various (V, P) — > 
(P, V)7 transitions, with P = rj, r]', taken from the PDG M and the corresponding predictions for 



gf P from Eqs. (M) —Fit 1— and (4.4) —Fit 2 



For these reasons, we have performed a new fit ignoring the two K* — » if 7 channels. This 
new fit obviously gives the same results for 5, and z q whereas m s /fh and z s always 
appear in the combination z s m/m s which is fitted to 0.63 ± 0.02, in agreement with our 
previous values in Eq. (|4.3|) . Therefore, our results are insensitive to eventual modifications 
of future and desirable new data on K* — * K'y transitions. 

In Table |l[ we present a comparison between experimental data for the relevant VP7 
transitions with P = rj,r]' and the corresponding theoretical predictions (in absolute value) 
calculated from the fitted values in Eq. Q4.3D - We do not include in that comparison 
neither V — > ttj nor K* — ► Kj modes since they are used to constrain the complementary 
parameters g , 4>v , m s /fh and zk- The agreement is very good and all the predictions 
coincide with the experimental values within la. Notice the small experimental error for 
the g$rn coupling as compared to the theoretical one, which, as we will see in Sec. ^, serves 
to highly constrain the allowed values for the r\-r{ mixing angle (pp. 

Now that we have performed a fit under the hypothesis of no gluonium we return to 
the main issue of this analysis, the phenomenological determination of the gluon content of 
the 77 and 7/ mesons. As stated before, a simultaneous fit of the overlapping parameters z qyS 
and the three mixing angles cj>p, 4> v g an d 4>n'G is n °t feasible. Therefore, we first assume 
<fir)G = 0) *- e - ^r) = 0, and then proceed to fit the gluonic content of the rj' wave function 
under this assumption. The results of the new fit are 3 

g = 0.72 ± 0.01 GeV" 1 , ^ = 1.24 ±0.07, 4> v = (3.2 ± 0.1)° , 

<j> P = (41.4 ± 1.3)° , \^ G \ = (12 ±13)° , (4.4) 
z q = 0.86 ± 0.03 , z s = 0.79 ± 0.05 , z K = 0.89 ± 0.03 , 

with x 2 /d.o.f.=4.2/4. The quality of the fit is similar to the one obtained assuming a 
vanishing gluonic admixture for both mesons (x 2 /d.o.f.=4.4/5). The fitted values for z q 



There is a sign ambiguity in 4> v 'g that cannot be decided since this angle enters into X n i and Y v > 
through a cosine. 
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and z s are compatible with those of Eq. (4.3). The result obtained for 4>wq suggests a very 



small amount of gluonium in the rj wave function, in fact compatible with zero within la. 
Using Eq. ( |A.6 ) to calculate Z^i from fi^c gives \Z v /\ = 0.2 ± 0.2. This is one of the main 



results of our analysis. Accepting the absence of gluonium for the r\ meson, the gluonic 
content of the rj wave function amounts to \^ G \ = (12 ± 13)° or Zj = 0.04 ± 0.09. 

In other words, our values for <f>p and (prj'G ( or ^n') contrast with those reported by 
KLOE recently <j) P = (39.7 ± 0.7)° and |c/y G | = (22 ± 3)° —or = 0.14 ± 0.04— 0. As 
indicated in Sec. |l], a possible explanation of this discrepancy could be the use in Ref. Q 
of old values for the overlapping parameters that the present analysis tries to update. In 
Table ffl, we also include the theoretical predictions for the various transitions involving r\ 



or rj calculated from the fitted values in Eq. (L4). As expected, there is no significant 
difference between the values obtained allowing for gluonium (Fit 2) or not (Fit 1) in the 
rj wave function. Likewise, we predict the value of the ratio 

_T(<f>-+rfi) 2 2 / m s z q tan <j> v \ 2 f p v > \ 3 

R <t> = r r(2 - = cot P cos <p v , G 1-— — . , — , (4.5) 

r(0^?77) V m z s sm2(j)p J \p v J 

to be (4.7 ±0.6) x 10~ 3 , in agreement with the experimental value in Ref. Q, (4.8 ±0.5) x 
10~ 3 , and the most recent measurement by KLOE g]], (4.77 ± 0.09 sta t ± 0.19 syst ) x 10~ 3 . 

For completeness, we perform another fit assuming from the beginning a null gluonic 
content for the rj meson. Consequently, we fix (ft^c = and leave free. The results 
obtained are the following: 

g = 0.72 ± 0.01 GeV- 1 , ^ = 1.24 ±0.07, <f> v = (3.2 ± 0.1)° , 

<Pp = (41.5 ± 1.3)° , |^g| - 0° , (4.6) 

z q = 0.86 ± 0.04 , z s = 0.78 ± 0.06 , z K = 0.89 ± 0.03 , 

with x 2 /d.o.f.=4.4/4. The fitted value for (j)„Q is very close to zero. For that reason, it 
is better to express this value in terms of the more common parameter. As a result, 
one gets Z^ = 0.00 ± 0.12, thus showing a vanishing gluonium contribution in the r\ wave 
function. This is a second important result of our analysis which complements the one 



discussed after Eq. (4.4). To sum up, the current experimental data on VPj transitions 
seem to indicate within our model a negligible gluonic content for the 77 and rj mesons. 

A final exercise we have done is to check whether the very recent measurements (not 
included in Ref. 0) on p,oj,(f) — > r}^ from the SND Coll. 0] and (f> — > r/7 from KLOE 
Q modify the results of our analysis. The values of the couplings associated to these new 
data are displayed in Table ||. As shown, the central values are nearly the same as those 
from Ref. Q, except for p — > r/7, whereas the errors for p,u>,(p ~^ are comparable to 
the world averages and the error for <j) — > r/7 is reduced by a factor of three. Assuming 
absence of gluonium in the r/ and r/ wave functions the results of the fit are 

(p P = (42.7 ±0.7)° , z q = 0.83 ± 0.03 , z s = 0.79 ± 0.05 , (4.7) 

with x 2 /d.o.f.=4.0/5. Accepting now a gluonic admixture only in the r( one obtains 

4> P = (42.6 ± 1.1)° , I (frn'G I = (5 ±21)° , z q = 0.83 ± 0.03 , z s = 0.79 ± 0.05 , (4.8) 
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Transition 


Sy£ 7 (latest) 


5^p 7 (Fit 3) 


4 h p 7 ( Fit 4 ) 


P ~ 


-> VI 


0.429 ± 0.023 


0.436 ±0.017 


0.437 ± 0.028 


rf- 


♦ A 


U.41 ± U.Uo ULi J 


n /in 4- n no 
U.4U ± u.uz 


n /in 4- n n/i 

U.4U ± U.U4 


UJ - 


■* VI 


0.136 ±0.007 


0.134 ±0.006 


0.134 ±0.009 


rf- 




0.139 ±0.015 (PDG) 


0.146 ±0.006 


0.146 ±0.013 


4>- 


-> r/7 


0.214 ±0.003 


0.214 ±0.017 


0.214 ±0.012 


<f>- 


* f/7 


0.216 ±0.005 


0.216 ±0.019 


0.216 ±0.018 



Table 2: The same as in Table | but for g$ p from Eqs. (g/7|) —Fit 3— and (|J) —Fit 4- 
compared to the latest <?yp 7 (in GeV -1 ) from Refs. [|, [T^ |. 



with x 2 /d.o.f.=4.0/4. The values of the remaining parameters are the same as in Eqs. 
and (fOJ), respectively. As one can see, the central values for the mixing angle <j>p slightly 
increase while the ones for c/) v 'g (or Z^, = 0.01 ± 0.07) and z q decrease. However, in both 
cases the quality of the fit is as good as it was when only the experimental data provided by 
Ref. 0] were considered. Moreover, the fitted values of all the parameters are compatible 



within la with those of Eqs. (4.3,4.4). The corresponding theoretical predictions for the 



gvP'y couplings with P = rj,rj' calculated from the fitted values in Eqs. (4/7) — Fit 3 — and 
Q4.g| ) — Fit 4 — are also shown in Table ||. Note in this case the small experimental error for 
g^n'-y which will constrain even more the allowed values for Z^, (see Sec. ||). In conclusion, 
the latest experimental data seem to confirm the null gluonic content of the r\ and rf wave 
functions. 



5. Comparison with other approaches 

Our main results can also be displayed graphically following Refs. Jl], ||, |(|. This will 
serve us to present the bounds obtained using our approach and compare them with other 
approaches. Very briefly, the analysis of Ref. || is based on the 577(3) quark model 
supplemented by the SU(3) breaking parameter fh/m s dealing with the difference of the 
down and strange quark magnetic moments but without taking into account the vector 
mixing angle <j)y or the different vector-pseudoscalar overlapping parameters, both specific 
features of our approach. The experimental analysis of Ref. [[l| is based on the same 
approach including the latter two features. Finally, the analysis of Ref. § utilizes the 
SU(3) quark model and the vector meson dominance model (VMD). In this framework, 
the transition amplitudes are expressed in terms of masses and decay constants affected by 
SU (3) breaking effects. 

In Fig. 0, we plot the regions for the X v and Y„ parameters which are allowed by the 
experimental couplings of the p — * 777, oj — > 777 and cj) — > 777 transitions (see Table ||). The 
limits of the bands are given at 68% CL or la. For the rest of the parameters involved in 
the determination of the allowed regions we have used the fitted values for g, (fiy, m s /fh, 
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Figure 1: Constraints on non-strange (X v ) and strange (Y v ) quarkonium mixing coefficients in 
the 77 as defined in Eq. ( |2.1| ). The circular boundary denotes the constraint X'^ + < 1. The 
mixing solutions corresponding to the 77 being a pure octet (X^ = — -t=Y" ?) = -^) — open circle — 
and (X v = = — closed circle — are shown. The straight lines for the upper and lower bounds 
of the mixing solution in absence of gluonium, <j)p = (41.5 ± 1.2)°, are also shown. The vertical 
and inclined bands are the regions for X v and Y v allowed by the experimental couplings of the 
(p, u>, 4>) — > 777 transitions in Table 111. 



z q , z s , and zk from Eq. (|4.6|) . In addition to the bands, we have also plotted the circular 
boundary denoting the constraint Xl + < 1 as well as the favoured region for the 
rj-rj' mixing angle assuming the absence of gluonium, 40.3° < (ftp < 42.7°, obtained at la 
from the corresponding fitted value in Eq. (4.3). As one can see, there exists a perfect 
intersection region of the three bands located precisely on the circumference at the region 
preferred by the mixing angle, thus indicating a vanishing gluonium contribution in the 77 
wave function, i.e. \Z V \ = 0. The small size of the intersection region, mainly due to the 
small experimental error in the coupling, and its precise location on the circumference 
gives strong evidence in favour of the former statement. Notice also the importance of 
considering the vector mixing angle (j)y different from zero, which translates into a finite 
slope for the (u, <p) — > 7/7 bands, in the result obtained. This contrasts with the analysis of 
Refs. H], |(§ where the vector mixing angle was not taken into account at that time because 
their effects were covered by the bigger experimental errors of the bands. With the present 
available experimental data, fixing the value of this angle to zero would have lead to an 
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Figure 2: Constraints on non-strange (X v i) and strange (Y n >) quarkonium mixing coefficients in 
the rj' . The mixing solutions corresponding to the 77' being a pure singlet (X v i = V^Y^i = 
— open circle — and {X,y — Y n > — -^) — closed circle — are shown. The vertical and inclined bands 
are the regions for X n i and Y v i allowed by the experimental couplings of the 7/ — > (p, uj)"f and 
4> ~ * v'l transitions. 

incompatible solution for the three bands. For completeness, the points corresponding 
to the 77 being a pure octet (X v = — -±=Y n = -^) and the "democratic" mixing solution 
(X^ = = 4|) are also shown in Fig. || as an open and closed circle, respectively. The 
latter solution, which is now excluded, was still acceptable in the analysis of Ref. ||, where 
the upper bound \Z V \ < 0.4 was obtained from the processes (p,(f>) — ► and 77 — ► 77. 
In Ref. after correcting some inconsistent data, there was no region where all three 
constraints, coming from the (u, <j)) — * 777 and 77 — > 77 decays, overlap. 

The situation is not so constrained for the 77' as we proceed to discuss. In Fig. ||, we 
plot the regions for the X„i and Y„i parameters which are allowed by the experimental 
couplings of the 77' — » pj, rf — > wy and <j) — > rf'y transitions. The remaining parameters 
are taken from Eq. ( |4.4| ). Again, we also plot the constraint implied by X"*, + Y3 < 1 
and the favoured region for the 77-77' mixing angle in absence of gluonium. This time there 
exists an intersection region of the three bands inside and on the circumference. As most 
of this region is interior but close to the circular boundary it may well indicate a small but 
non necessarily zero gluonic content of the 77'. Indeed, we have found Zl = 0.04 ± 0.09 
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Figure 3: The gray ellipse in the (</>p, Z^,) plane corresponds to the allowed region at 68% CL of the 
solution in Eq. (4.4) — left plot — and Eq. ( |4.S| ) — right plot — , respectively, assuming the presence of 
gluonium. The different bands are the regions for tfip and Z^, allowed by the experimental couplings 
of the 77' — > p7 (dashed line), 77' — > 077 (dot-dashed line), <ri — * 777 (dotted line), and — > 77' 7 (solid 
line) transitions in Table [l| — left plot — and Table ^| — right plot — , respectively. 



(or \Z V '\ = 0.2 ± 0.2) or using the angular description \4>j 1 'g\ = (12 ± 13)°. The size of 
the error is precisely what prevent us from drawing a definite conclusion concerning the 
amount of gluonium in the 77' wave function. More refined experimental data, particularly 
for the <p — ► 77' 7 channel, will contribute decisively to clarify this issue (see below). Clearly, 
the inclusion of this process is of major importance for the determination of the gluonic 
admixture in the 77', as observed for the first time in Ref. g. In this latter analysis, 
where 77' — > pj and 77' — > 77 were used in addition to different values for cj> — ► 7/7, the 
absence of a significant constraint on Yw was keenly felt. However, the mixing solution 
Xrf =Y V ' = -^ — closed circle in Fig. ^ — was still acceptable. It was not the case for the 
solution identifying the r( as a pure singlet, = y/2Y v i = 4-, — open circle in Fig. || — . In 
the present analysis, the "democratic" mixing solution is excluded at the la level whereas 
the singlet solution is clearly excluded. In Ref. Q, where 77' — ► 0*7 was also included in the 
analysis, the maximum gluonic admixture in the 77' was obtained to be 26% for 6p = —11° 
(or <j)p ~ 44°). In our case, for = 0.04 ± 0.09, one gets R = x ,+j' l+z ■ = (13 ± 13)% 
for <j)p = 41.4° (or Op = —13.4°). The same analysis also anticipated that the existence of a 
gluonic content for the 77' would be excluded for large \0p\ (or small 4>p). Finally, the KLOE 
analysis in Ref. Q has found a solution allowing for gluonium with aV = (39.7 ±0.7)° and 
Z^i = 0.14 ±0.04 — or |<^yc?| = (22 ±3)° — , in disagreement, particularly for the Zy value, 
with our results. A second analysis without the constraint from the 77' — ► pj decay, which 
is the less precise, gives (ftp = (39.8 ± 0.8)° and Z£ = 0.13 ± 0.04, thus not modifying 
their conclusions. To make this difference more graphical, we follow Ref. 0] and plot in 
Fig. H — left plot — the constraints from if — > (p,uj)j and (f> — > (rj,rf)~f in the (4>p,Z^,) 
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plane together with the 68% CL allowed region for gluonium as obtained from Eq. (|4,4|), 
The point corresponding to the preferred solution, (cj)p,Z^,) = (41.4°, 0.04), is also shown. 
The allowed region is very constrained in the (ftp axis by the experimental value of the 
g^ rn coupling, whose vertical band denotes its non dependence on Z^,. The other three 
bands, all dependent on (ftp and Z^,, constrain the amount of gluonium down to a value 
compatible with zero within la. As mentioned before, our result differs from the KLOE's 
one, where (</>p,Z^,) = (39.7°, 0.14) is the preferred solution and the allowed region for 
gluonium is far from the <f>p axis (see Fig. 5 in Ref. Q). The value compatible with zero 
we have obtained for the gluonic content of the rf wave function is fully confirmed from a 
graphical point of view as soon as one includes the latest measurements of the eft — > (rj, 7/) 7 
decays from Refs. ]l], E3. As seen in Fig. ||| — right plot — , the allowed region for gluonium 
decreases due to the smaller experimental error of the g^'-y coupling (although the error 
of g^y increases a little) and its central point (the preferred solution) is now located at 
(cj)p,Z^,) = (42.6°, 0.01), even much closer to the (j)p axis than before. 

In all the other approaches discussed so far, the decays (r/, rf) — > 77 were also included 
in the analyses. As shortly stated in Sec. |l[ in our approach we do not take into account 
these processes since they are well explained using a two mixing angle scenario M, M. 
Nevertheless, our framework can be extended to describe P° — > 77 decays in a similar way 
to V — > P7 and P — > V7. In such case, the needed ratios of decay widths are found to be 



J^ = Sfe) 8 ( 5 ^ Jf " + ^ 5 - y ") 2 



where, in addition to the parameter rh/m s related to the different quark spin-flip effects 
of the emission of one photon, we have to introduce two new SU(3) breaking parameters, 
z q = C q /C n and z s = C s /C n , dealing with the spatial wave functions effects of quark- 
antiquark annihilation into the other photon. These new parameters are analogous to 
the overlapping parameters, z q and z s , of the VP7 transitions, and for the same reason, 
the U(1)a anomaly, we also distinguish between the quark- antiquark annihilation effects 
in |vr) and \rj q ). For comparison, in Ref. H the former SU(3) breaking effects are not 
considered, i.e. z q = z s = 1 and fh/m s = 1, while in Ref. ||, where the P° — > 77 
decays are characterized by means of pseudoscalar decay constants, one identifies = z7. x 
and t 2 - = ^-zT 1 - As seen in Eq. (|5.1|), these ratios are not useful for fixing, within 

J TV ill " * 

our approach, the mixing parameters of the rj-rj' system since they also depend on the 
unknown values of z q and z s . However, in order to check the consistency of this extended 
framework, one could use for the mixing parameters and fh/m s the results obtained from 
the analysis of V — » P7 and P — > V7 decays and then fix the values of z QyS from the two 
ratios under consideration. In this case, one gets z q = 0.96 ± 0.03 and z s = 0.79 ± 0.13 
and hence f q = (1.05 ± 0.03)/^ and f s = (1.57 ± 0.28)/^, in agreement with the values 
f q = (l.lOiO.03)/,,- and f s = (1.66±0.07)/,r found in Ref. j?]], thus showing the consistency 
of the mixing parameters obtained and of the whole approach. 



- 12 - 



6. Summary and conclusions 

In this work we have performed a phenomenological analysis of radiative V — > -P7 and 
P — > V"f decays with the purpose of determining the gluon content of the 77 and rj mesons. 
The present approach is based on a conventional SU(3) quark model supplemented with two 
sources of SU(3) breaking, the use of constituent quark masses with m s > fh, thus making 
the VP"f magnetic dipole transitions to distinguish between photon emission from strange 
or non-strange quarks, and the different spatial extensions for each P or V isomultiplet 
which induce different overlaps between the P and V wave functions. The use of these 
different overlapping parameters — a specific feature of our analysis — is shown to be of 
primary importance in order to reach a good agreement. 

Our conclusions are the following. First, the current experimental data on VPj tran- 
sitions indicate within our model a negligible gluonic content for the r] and rf mesons, 
Z^ = 0.00 ± 0.12 and Z^, = 0.04 ± 0.09. Second, accepting the absence of gluonium for 
the 77 meson, the gluonic content of the 7/ wave function amounts to \4> v 'g\ = (12 ± 13)° 
and the 77-7/ mixing angle is found to be (f)p = (41.4 db 1.3)°. Third, imposing the absence 
of gluonium for both mesons one finds <f>p = (41.5 ± 1.2)°, in agreement with the former 
result. Fourth, the latest experimental data on (p, to, <fi) — ► 777 and <fi — ► 7/7 decays confirm 
the null gluonic content of the r] and 77' wave functions. Finally, we would like to stress 
that more refined experimental data, particularly for the 4> v'l channel, will contribute 
decisively to clarify this issue. 
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A. Euler angles 

In presence of gluonium, the wave functions of the 77 and 77' mesons and the glueball-like 
state 1 can be decomposed as 

\ V ) = X v \ Vq ) + Y v \ Vs ) + Z v \G) , 

| ? /) = X v ,\ Vq ) + Y v ,\ Vs ) + Z V ,\G) , (A.l) 
\l) = X^) +Y L \r ]s ) + Z L \G) , 

where |?7 g ) = ^\uu + dd), \r/ s ) = \ss) and \G) = |gluonium). The t or 7/(1440) state, 
which we refrain from discussing here, could be identified with the 7/(1405) pseudoscalar 
resonance (see Ref. [|j for details). The nine coefficients in Eq. ( |A.1| ) are constrained by 
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three normalization conditions 



X 2 + y, 2 , + ZJ = 1 , (A.2) 



and three orthogonality conditions 

XjjXrf + YrjYrf + Z^Z^i = , 

x„x t + y„y t + z v z L = o , (A.3) 

X v ,X b + Y v ,Y L + Z v ,Z, = . 

Altogether implies that only three independent parameters are required to describe the 
rotation between the physical states (77, rf and t) and the orthonormal mathematical states 
(r] q , T] s and G). The rotation matrix can be written in terms of three mixing angles, <pp, 
4>rjG and 4> v >g-i which would correspond to the three Euler angles for a rotation in real, 
three-dimensional space. Explicitly, 

S^pC^^iQ - C^pS^icS^nG cfypCfyrfG + S^pS^^iQS^nG -s<f> v i G c<f> v G 1 I Vs I ) (A. 4) 
S<f) P S(f> v iG + C<t>pC<t> v 'GS(firiG C^pS^iQ - S<f> P C<f> v > G s 4>vG C<f> v >G C( PvG J \G J 

with (c, s) = (cos, sin). In the limit in which fi^a = fyrfG = 0, the gluonium decouples and 
the former matrix reduces to the usual rotation matrix describing 77-77' mixing, 





cos <pp — sm<pp 
sin0p cos (f)p 




(A.5) 



where (f)p is the 77-77' mixing angle in the quark- flavour basis related to its octet-singlet basis 
analog through Op = <fip — arctan\/2 ~ (f>p — 54.7°. An interesting situation occurs when 
the gluonium content of the 77 meson is assumed to vanish, i.e. (fi v c = 0. In this particular 

case, 

X„ = cos <pp , Y„ = — sin^p , Z„ = , 

(A.6) 

Xrf = sm(f)p cos 4>n'G > = cos <Ap cos <?Vg > ^7?' = — sm< ?VG • 

References 

[1] F. Ambrosino et al. [KLOE Collaboration], arXiv:hcp-ex/0612029. 

[2] A. Aloisio et al. [KLOE Collaboration], Phys. Lett. B 541 (2002) 45 [arXiv:hep-ex/0206010]. 

[3] A. Bramon, R. Escribano and M. D. Scadron, Phys. Lett. B 503 (2001) 271 
[arXiv:hcp-ph/0012049]. 

[4] W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 

[5] J. L. Rosner, Phys. Rev. D 27 (1983) 1101. 

[6] E. Kou, Phys. Rev. D 63 (2001) 054027 [arXiv:hep-ph/9908214]. 



- 14 - 



[7] R. Escribano and J. M. Frerc, JHEP 0506 (2005) 029 [arXiv:hep-ph/0501072]. 

[8] T. Fcldmann, Int. J. Mod. Phys. A 15 (2000) 159 [arXiv:hep-ph/9907491]. 

[9] G. Li, Q. Zhao and C. H. Chang, arXiv:hcp-ph/0701020. 

[10] N. Isgur, Phys. Rev. Lett. 36 (1976) 1262. 

[11] P. J. O'Donncll, Rev. Mod. Phys. 53 (1981) 673. 

[12] S. I. Dolinsky et al, Phys. Rcpt. 202 (1991) 99. 

[13] A. Bramon, R. Escribano and M. D. Scadron, Eur. Phys. J. C 7 (1999) 271 
[arXiv:hcp-ph/9711229]. 

[14] M. N. Achasov et al, Phys. Rev. D 74 (2006) 014016 [arXiv:hep-ex/0605109]. 



- 15 - 



